Andreev reflection tuned by magnetic barriers in superconductor contacted graphene 
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The Andreev reflection of the normal state-superconductor junction both in monolayer and bi- 
layer graphene with a single magnetic barrier is investigated by means of the Green's function 
formalism. Within the tight-binding model, we study the direction-dependent Andreev reflection 
of two-dimensional graphene-superconductor junctions in the specular and retro-reflection regimes. 
The presence of a magnetic barrier close to the superconducting hybrid junction introduces a rich 
phenomenology. Such a barrier is capable of tuning the preferred angles of incidence for the An- 
dreev reflection. In particular, it can enhance the specular reflection probability for certain angles 
of incidence in bilayer-based hybrid structures. When transmission is permitted, the Andreev re- 
flection manifests itself in isolated peaks and asymmetric resonances associated with offsets and 
Fano-type oscillations in the transmission, respectively. Moreover, Fabry-Perot oscillations in the 
Andreev reflection due to the interior scattering inside the magnetic barrier may appear. The im- 
pacts of magnetic barriers on the monolayer and bilayer hybrid interfaces are furthermore studied 
by calculating the differential conductances within the Blonder-Tinkham-Klapwijk formula. 

PACS numbers: 74.45.+C, 72.80.Vp, 73.23.-b 
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I. INTRODUCTION 

A dissipative electric current is converted into a dissi- 
pationless supercurrent by the celebrated Andreev retro- 
reflection at the normal-superconducting (N-SC) inter- 
face where an electron is reflected into a hole which re- 
traces the path of the incident electron^ In graphene, it 
has been suggested that due to its peculiar band struc- 
ture, the Andreev reflection (AR) can be specular for 
excitation energies above the Fermi energy. The specu- 
lar AR in graphene differs from the usual retro-reflection 
with clear experimental signatures^ Contrary to two- 
dimensional graphene, the AR in graphene nanorib- 
bons (GNR) is sensitive to the ribbon width and the 
pseudoparity of quantum states^ Besides the monolayer 
graphene (MLG)-superconductor junction, the bilayer 
graphene (BLG)-supcrconductor junction has been stud- 
ied by Ludwig, and a remarkable suppression of the dif- 
ferential conductance at voltages just below the super- 
conducting gap was founds 

Beyond a simple N-SC junction, a number of extended 
hybrid structures have also been proposed and studied, 
like for example the graphene-insulator-SCS and the SC- 
graphene-SC systems^ as well as a multi-terminal system 
of graphcnc-SC junctions where the subgap structure and 
crossed Andreev reflection were discussed^. In partic- 
ular, the transport properties through the graphene-SC 
junctions in the presence of magnetic fields have been 
studied. The distinct transport regimes are recognizable 
in the phase diagram^ 

On the other hand, the intriguing properties of 
graphene with normal contacts have been extensively 
investigated both experimentally and theoretically.— ~— 
Graphene exhibits a variety of exotic transport behav- 
iors, such as the anomalous integer and fractional quan- 
tum Hall effect and minimal conductivit y 1 1 ' 16 , as well 
as, an extremely high mobility due to the fact that the 



Dirac fermions can pass perfectly through classically for- 
bidden region (Klein tunneling )^2r—. This effect, how- 
ever, impedes the the confinement of Dirac fermions for 
the purpose of future applications. It has been shown 
that a local inhomogeneous magnetic field, i.e. a mag- 
netic barrier (MB), is able to restrict massless Dirac elec- 
trons effectivelyj 21 ' 22 Moreover, the electronic tunneling 
through various magnetic barrier geometries in graphene 
has been studied systematically by Masir et a h 23 ' 24 These 
studies have paved the way for constructing direction- 
dependent wave vector filters based on the MB which is 
viewed as a key element for graphene-based electronics^ 
In view of the fact that superconducting contacts as well 
as MBs in graphene give rise to numerous novel trans- 
port phenomena, it appears to be a natural next step to 
combine these two components and explore the resulting 
transport properties. In terms of fundamental physics, 
such a study will deepen our insight into the impact of 
MBs on the interplay between superconductivity and rel- 
ativistic dynamics. Furthermore, it will also provide a 
guideline for possible future experimental steps in this 
direction. 

In this paper, we investigate the AR both at the MLG- 
SC and the BLG-SC junctions with a single MB by utiliz- 
ing the real-space Green's function formalism. The angu- 
lar dependencies of the AR at the electron energy lower 
or slightly higher than the superconducting energy gap 
are calculated. The response of the hybrid junction for 
a single MB is examined in different transport regimes. 
Also, the differential conductance of the NSC interface 
with a MB closeby is considered to facilitate comparison 
to experiments. 

This paper is organized as follows. In Sec. II wc sketch 
the considered geometry and introduce the formalism. 
This is followed by the presentation and discussion of 
the numerical results in Sec. III. The summary and con- 
clusions constitute Sec. IV. 
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FIG. 1. (Color online) Schematic geometry of the hy- 
brid structure under consideration in monolayer or bilayer 
graphene. The central graphene device, where a local in- 
homogeneous magnetic field is applied, is connected by two 
semi-infinite graphene leads. The right lead is contacted by 
a superconducting metal to induce superconductivity. The 
tight-binding model satisfies boundary conditions to model 
two-dimensional graphene. Lower panel: A typical vector po- 
tential profile for a single magnetic barrier. 



II. MODEL AND FORMALISM 

We consider the two-terminal hybrid graphene- SC ge- 
ometry sketched in Fig. [I] The left and right leads 
are made of the semi-infinite metallic and superconduct- 
ing graphene, respectively. In the central region, a sin- 
gle magnetic barrier is applied to the normal-conducting 
graphene. This system can be described by the Dirac- 
Bogoliubov-dc Gennes (DBdG) Hamiltonian with finite 
and vanishing pair potentials A = Aoe l</3 corresponding 
to the superconducting and normal-states graphene, re- 
spectively. In this work we consider the s-wave symmet- 
ric pair potential, i.e., Ao is real and the superconducting 
phase tp = 0. 

The easiest way to induce superconductivity in 
graphene relies on the proximity effect by superimpos- 
ing a superconductor (SC) on graphene as shown in the 
upper panel of Fig. [1] Recent experiments have observed 
unambiguously superconductivity in graphene via this 
mechanism<2£r— The details of the induction of super- 
conductivity and the lattice matching between graphene 
and the SC are of no interest for our purpose. 

We make use of the tight-binding (TB) approach for 
the description of the transport calculations. The DBdG 
Hamiltonian for MLG in TB model reads 
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where c\{ci) is the creation (annihilation) operator on 
the site i. The first sum runs over all the atomic sites 
with on-site electrostatic potential Vi. Ep is the Fermi 
energy, and the 2x2 Pauli matrices a = (a x , a v , a z ) 
act on the electron-hole space. a± — a x ± ia y are the 



raising and lowering operators, respectively. The pair 
potential A(x) describes the coupling between electrons 
and holes in superconducting states, which vanishes in 
the graphene normal state. The second sum runs over 
all the pairs of neighboring atomic site with the 

hopping energy 70 ~ 3 eV. All states are spin-degenerate. 

The Hamiltonian describing the BLG in Bernal stack- 
ing is written a; 
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where %mlg * s ^ nc Hamiltonian of the v-th graphene 
layer given by Eq. ([T]), c vi (d Vt i) is the creation (annihi- 
lation) operator at sublattice A(B), in the layer v = 1, 2. 
71 = 0.39 eV is the coupling energy between sublattice B 
and A in different graphene layers, which is the most rel- 
evant for the low-energy band of BLG. Other inter layer 
coupling energies, in particular 73 and 74 leave the basic 
phenomenology unaltered and are therefore not consid- 
ered here i 31 i 32 

The transport properties of the two-terminal system 
is calculated by means of the real-space Green's func- 
tion formalism which has been discussed in detail in Refs. 



[33ll34| . Here, we merely give the most relevant formalism 
for the present calculations. The total retarded Green's 
function is calculated by the recursive algorithm and con- 
nected to the two leads by 



G R (e) = (E F + e-Hi 



(3) 



where e is the excitation energy relative to the Fermi en- 
ergy Ep. T-Ld is the Hamiltonian of the central device 
region including scattering, the self-energies ^l(r) a c- 
counts for the effects of the semi-infinite left and right 
leads, respectively. The right lead is superconducting, 
while the left lead where the- wave enters is always in the 
normal state. The transmission amplitude for a particle 
p in mode n entering from the left lead and exiting into 
the right lead as the particle p' in mode m is given by 
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where G% +h0 = (N + l\G R (e)\0) is the Green's function 
clement in the real-space with and N + 1 denoting the 
first slice in the left and right leads, respectively, IVm 

is the spectral matrix given by ^l(r) = i(^i(R) — ^;(m)- 
0n(m) is the eigenfunction of the Bloch state n(m) with 
the group velocity v n(m) . 

The reflection amplitude of a particle p in mode n in 
the left lead I to a particle p' which may be either an 
electron (e) or a hole (h) of the mode m in the same lead 
I reads 
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where the Green's function element Gq — (0\G R (e)\0). 
Apparently, r he and r ee denote the Andreev and normal 
reflection amplitudes, respectively. 

The corresponding reflection R and transmission T 
probabilities are given by 



R(T) 
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Henceforth we use R A and R N to denote the Andreev 
and normal reflections, respectively. 

By imposing periodic boundary conditions on the 
transverse (y-) direction (thereby wrapping the chain 
into a loop and removing its ends), the Green's func- 
tion technique is able to treat the transport properties 
of 2D graphene layers. In our calculations, this is imple- 
mented by connecting the upmost and downmost carbon 
atoms of the GNR of a width of W by the hopping energy 
j e lkyW with the wave vector k y in the y-direction. This 
technique is rather efficient for dealing with the transport 
properties of the 2D graphene^ 

To incorporate the perpendicular magnetic field in 
this model, we start with the Dirac equation instead 
of Peierls' substitution usually applied within the TB 
model. In more detail, the momentum in the presence of 
magnetic fields becomes p + eA(x). By using the Landau 
gauge, the vector potential is given by A = (0, B(x)x, 0) 
with B(x) being the inhomogeneous magnetic field along 
the a;-direction. In our calculations, the MB is parame- 
terized by its length L and strength B, i.e. a rectangular 
barrier, which generates a vector potential profile shown 
in the lower panel of Fig. [1] A more realistic barrier 
shape as implemented in experiments^ is shown in the 
main figure, but the shape modifies the vector potential 
profile only slightly in its incremental part and does not 
affect our numerical results significantly. Consequently, 
the MB is incorporated by rewriting the periodic hopping 
energy as ^ e ^k v +A y {x)}w w ^ A y (x) = B(x)dx. 

The differential conductance of the system is calculated 
according to the Blondcr-Tinkham-Klapwijk formula^ 
as 



9 = 



QT ft/ 2 

^77 = .90 / (1 - R N + R A ) cos ada, (7) 

OV J-tt/2 



where go is the ballistic differential conductance for cither 
MLG or BLG and includes the spin and valley degenera- 
cies. The normal reflection R N and Andreev reflection 
R A are calculated from Eq. (JBJ. 



III. RESULTS AND DISCUSSION 

In this section we present and discuss the numerical re- 
sults of our transport calculations. The hybrid structures 
in both MLG and BLG are investigated. We consider two 
different regimes, i.e. the subgap regime with incoming 
energy e < Aq and the open regime with e > Aq. 



A. subgap regime e < Ao 

We first focus on the subgap regime. The AR in this 
regime can be further identified as the retro- (e < Ep) 
and specular (e > Ep) reflection depending upon the 
excitation energy relative to the Fermi energy. Fig. [2] 
shows the AR as a function of the angles of incidence with 
increasing of the strength of MBs. The AR vanishes if 
the angle of incidence is beyond a critical angle a c which 
is given by 



E f 



(8) 



with v — 1,2 for monolayer and bilayer graphene, 
respectively*^ Inspection from the dispersion relation, 
this is because a gap is opened in the hole band for small 
longitudinal momentums. The numerical results coincide 
well with the Eq. In the absence of MBs, a com- 

mon feature for the retro-reflection and specular reflec- 
tion in MLG is the unity reflection at normal incidence 
(a = 0) due to the conservation of chirality in contrast to 
the usual NSC junctions. Namely, the reflection from an 
electron to a hole at the graphcne-SC interface does not 
involve the conversion of sublattices, whereas the normal 
reflection R N does involve scattering between sublattices 
at a = 0£ Note that the AR probability of 2 in our case 
is due to the valley degeneracy. 

The magnetic field filters the wave vectors prior to the 
electron approaching the the NSC interface, making the 
AR probabilities asymmetric with respect to the normal 
incidence owing to the time-reversal symmetry breaking. 
As the strength of the MB increases, the suppression of 
AR by MBs becomes more prominent. For MLG hybrid 
structures, the retro- and specular reflections behave sim- 
ilarly. 

In comparison to the case of MLG, the retro-reflection 
and specular reflection in the BLG system exhibit much 
more distinct features as shown in Fig. [3] Considering 
that the dispersion of BLG is not perfectly parabolic, 
the numerical values of a c deviate slightly from those 
from Eq. ([8]). More significantly, the retro-reflections 
at e < Ep always take place with a certain probabil- 
ity at normal incidence, while the specular reflections at 
e > Ep always show zero probabilities corresponding to 
perfect normal reflections, i.e. R N = 1. This is consis- 
tent with the theoretical results based on Dirac theory 
and can be interpreted by the compatibility of the spinor 
wave functions of excitations in graphene^ To be more 
specific, the wave function of an incident electron in the 
sublattice basis is orthogonal to that of the reflected hole 
in the valence band so that specular AR is forbidden at 
a = 0. This feature persists in the presence of MBs, while 
the angular dependence of the AR is bent downwards. 

Furthermore, the MBs may raise the AR probability 
up to unity for certain directions, particularly for spec- 
ular reflections. The two wings of the specular AR are 
bent and expanded for the smallest magnetic field. As 
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(a) retroreflection (b) specular reflection 
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FIG. 2. (Color online) The angular dependence of the An- 
dreev reflection for monolayer graphene-SC hybrid junction 
with various MBs in the subgap regime, (a) The retro- 
reflection with e/E F = 0.05 and E F /A = 10, a c « 65°. The 
blue, red, and green curves from exterior to interior respec- 
tively correspond to Mb = 1.3, 0.9, and 1.3. (b) The specular 
reflection with e/E F = 19 and E F /A = 0.05, a c « 64°. 
The blue, red, and green curves from exterior to interior re- 
spectively correspond to Mb = 1.7, 1.2, and 0.9. For MLG, 
hv F k = E F + e with v F the Fermi velocity. 
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FIG. 4. (Color online) The Andreev reflection and corre- 
sponding transmission probabilities as a function of angles 
of incidence a in MLG with various MBs at e > Ao. (a) The 
retro-reflection regime with e/E F — 0.11 and E F /Ao = 10 
with the corresponding transmission (b). (c) The specular re- 
flection regime with e/E F = 5.1, and E F /Ao — 0.2 for various 
d with the corresponding transmission (d). The parameters 
for the MB in this case are fixed as L/Ib = 0.5 and Mb = 0.7. 



(a) retro-reflection 



(b) specular reflection 




cases. However, in other three cases (retro- and specular 
AR in MLG, and retro- AR in BLG), the AR probabilities 
nearly saturate for most of incoming angles such that the 
increments in ARs arc negligible. In the time-reversed 
N-SC junctions, elastic scattering at the hybrid inter- 
face and the random impurities produces a constructive 
interference between the incident electron and the retro- 
reflected hole resulting in an enhancement in the differen- 
tial conductance— ~— This is known as the "reflectionless 
tunneling" and possibly responsible for the zero-bias con- 
ductance peak in graphenc hybrid junction— However, 
it is not a major effect here since the specular reflection 
and the magnetic field spoil the phase conjugation. 



FIG. 3. (Color online) The angular dependence of the An- 
dreev reflection for bilayer graphene-SC hybrid junction at 
various MBs in the subgap regime, (a) The retro-reflection 
with e/E F = 0.03 and E F /A = 2, a c « 70°. The blue, read, 
and green curves from exterior to interior respectively corre- 
spond to Mb = 2, 2, and 1.5. (b) The specular reflection with 
e/E F = 37 and E F /A = 0.02, a c » 71°. The blue, red, and 
green curves from exterior to interior respectively correspond 
to M B = 1-8, 1.8, and 1.3. For BLG, h 2 v 2 F k 2 / 7l = E F + e 
with vf the Fermi velocity. 

the magnetic field increases, the lower wing shrinks, while 
the upper wing persists. The enhancement in the upper 
wing of the AR can be understood as follows. The elec- 
trons may undergo successive scattering inside a resonant 
cavity formed by the MB and the N-SC interface. Each 
of the cycles can build up the probability of conversion 
from electrons to holes. This mechanism and the filter- 
ing effect of MBs coexist and compete each other in all 



B. Open regime e > Aq 

We now consider the case of incident electrons with 
an energy e slightly above the superconducting gap. In 
this case, normal transmission emerges and the AR pro- 
cesses are severely suppressed. In Fig. 2] we show the 
Andreev reflection with different MBs and d as a func- 
tion of angles of incidence a. The ARs without MBs are 
also plotted (the dashed curve) for comparison. Stronger 
MBs confine the AR successively towards a more nega- 
tive and narrower regime of angles of incidence because 
of the filtering effects of the MBs. In Fig. H (a), the 
distance d between the MB and the SC interface is also 
varied. An AR peak appears on the right hand side of the 
curve and diminishes as d increases. The local AR min- 
imum nearby the peak corresponds to the transmission 
"peak" through the N-SC junction as shown in Fig. HJb), 
which is a remnant of the protrusion in the transmission 
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FIG. 5. (Color online) The Andreev reflection probabilities 
as a function of angles of incidence a in BLG with vari- 
ous single MBs at e > Ao. (a) The retro-reflection regime 
with s/Ef = 0.11 and Ef /Ao = 10 with the correspond- 
ing transmission (b). (b) The specular reflection regime with 
e/Ef = 4.2, and Ef /Ao = 0.25 for various d with the cor- 
responding transmission (d). The parameters for the MB in 
this case are fixed as L/Ib = 1.0 and Mb = 2.3. 



without MBs (dashed curve). 

Fig. @|c) shows the AR probability in the specular 
regime with e > Ep, where we focus on the influence of 
d on the AR with a constant MB. Compared with the 
retro-reflection, a concave specular AR shows up around 
a = 0, accompanied by a convex-shaped transmission 
as indicated by the dashed curve in Fig. EJd). In the 
presence of MBs, the overall AR is suppressed except a 
small angular interval for the largest d. Interestingly, 
there is a clear crossover around a = —0.1 (in units of 
7r/2) separating the specular AR into two regions. In the 
right region (a > —0.1), the transmissions arc completely 
suppressed for all MBs. A possible scenario here is that 
the electrons are multiply scattered between the MB and 
the N-SC interface (see the inset of Fig. 0Jd) for an 
example), which compensates the loss in the AR due to 
the prior scattering of the MB to some extent. Therefore, 
the AR probability drops as the separation d between 
the MB and the superconducting lead increases. For the 
angles of incidence a < —0.1, it crosses over to another 
region, where the transmission emerges and Fano-type 
resonances can be observed as shown in Fig. |Ud). What 
happens here is that, the electron acquires an additional 
phase after the multiply scattering between the MB and 
the SC, then transmits into the superconducting lead. 
The additional phase results in aymmctric line shapes in 
terms of interference phenomenal Correspondingly, the 
resulting reflection resonances can be recognized at the 
angle a = -0.47. 

In Fig. ODJa), it is shown that the behavior of retro- AR 
of BLG-SC junctions is very similar to that of the MLG- 
SC junction. For d = 51b, the reflected holes may un- 
dergo interior multiple scattering inside the MB, thereby 



leading to the Fabry-Perot-type oscillations which also 
appear in the transmission of GNRs subject to a sin- 
gle MB and carbon nanotubes with normal contactsi 34 i 43 
For the specular reflection, the MB shifts the reflection 
minima to negative angles of incdience and raises the re- 
flection probability with the angles of incidence a > 0.2 
markedly as in the subgap regime. The growth of AR 
gets less prominent, arising from suppression of succes- 
sive reflections both from the MB and the SC interface 
as d increases. For a < 0.2, a similar crossover to an- 
other region with nonzero normal transmissions can be 
identified. 



C. Differential conductance 

We now turn to a discussion of the differential conduc- 
tance g of N-SC junctions. Fig. [6^a) shows g with various 
Ep I Ao for MLG hybrid junctions with and without MBs. 
The conductance drops to zero at eV = Aq since the hole 
states are unavailable at this voltage. The similar feature 
can also be found in the case of BLG. The effect of MBs 
exhibits by reducing the overall differential conductance 
for all applied voltages except at the singular points, al- 
beit the angular dependence of the AR has been washed 
out by the angular integration. 

The results for the BLG case are shown in Fig. [D^b). 
By applying the identical MBs on the BLG-SC junction 
for different Ep/A , we find that the MBs in BLG hy- 
brid structures are more transparent than those in MLG. 
Particularly, for Ep/A = 0.2, the conductance with a 
MB can be slightly higher than that without a MB, be- 
cause the MB enhanced the Andreev reflections primar- 
ily in the specular regime. For Ep/Aq > 1, the specular 
regimes are already outside the subgap regime and enter 
the region where the transmissions are dominant. As a 
result, the effect is less pronounced. Moreover, a conduc- 
tance dip, which arises from the zero AR at the normal 
incidence in the specular regime, shows up at the voltage 
just below the superconducting energy gap. This con- 
ductance dip is also found in the analytical result based 
on the Dirac equation^. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have studied the Andreev reflection at 
normal-conducting-superconduction graphene interfaces 
with magnetic barriers closeby. The magnetic barrier 
constrains the permissible angles of incidence, thereby 
leaving the Andreev reflection asymmetric with respect 
to the normal incidence, in contrast to electrostatic bar- 
riers. In the subgap regime, the magnetic barrier in bi- 
layer graphene shows a quite different behavior for the 
retro- and specular Andreev reflection compared with the 
monolayer hybrid junction. In particular, the magnetic 
barrier can enhance the Andreev reflection for specific 
angles of incidence due to the multiply reflection of clcc- 
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FIG. 6. (Color online) The differential conductance with a 
single MB as a function of applied bias for various E F /Ao 
ratios, (a) k F l B = 0.24,0.12,1.2 and L/l B = 0.5,0.2,0.5 
for MLG. Here hv F k F = E F . (b) k F l B = 2,1.4,4.5 and 
L/l B = 0.5,0.2,0.5 for BLG. Here t^v^kp/^ = E F . The 
dashed curves correspond to the cases without MBs. 

trons between the magnetic barrier and normal state - 
superconducting interface. 

In the open regime where the transmission is dominant, 
these hybrid junctions show a rich transport phenomenol- 
ogy. For the retro-reflection, the barrier deflects the prop- 



agation of the electron such that the transmission protru- 
sion penetrates into the Andrccv reflection regime, which 
is manifested by Andreev reflection "peaks" . The multi- 
ple scattering of electrons inside the barrier may be signif- 
icant at small angles of incidence and lead to Fabry-Perot 
interference. In the case of specular reflection, the An- 
drccv reflection can be identified into different transport 
regions, in which Fano-type resonances may appear. 

The magnetic barrier strongly suppresses the differen- 
tial conductance of monolayer hybrid structures, whereas 
it shows rather transparency in bilayer. One reason is 
the nearly parabolic dispersion of bilayer. In graphene 
contaced by a normal metal, electrons with a parabolic 
spectrum have shown robustness to magnetic barriers 
compared to those with the linear spectrum. Also, the 
magnetic barrier further intensifies the conversion from 
electrons to holes at the bilayer hybrid interface in the 
specular regime. In addition, as a result from the An- 
dreev reflection minimal at the normal incidence in the 
specular regime, the differential conductance of the bi- 
layer junction shows a dip just below the energy gap, 
which preserves in the presence of magnetic barriers. 

The phenomenology of the Andreev retro-reflection 
can be viewed as an electronic analogue of optical phase 
conjugation^ The presence of the magnetic barrier 
strongly modifies the preferable direction of the incident 
electron, and offers a possibility to build a wave-vector- 
dependent phase-conjugating device. 
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